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Abstract 

We compute corrections of relative order v'^ to the rates for the decays of ^Sq heavy quarkonium 
into two photons and into hght hadrons and for the decays of "^Si heavy quarkonium into a lepton 
pair and into hght hadrons. In particular, we compute the coefficients of the decay operators that 
have the same quantum numbers as the heavy quarkonium. We also confirm previous calculations 
of the order- 1;^ corrections to these rates. We find that the v expansion converges well for the decays 
of ^Sq heavy quarkonium and for the decay of heavy quarkonium into a lepton pair. Large 
higher-order-in-t; corrections appear in the decay of heavy quarkonium into light hadrons. 
However, we find that the series of coefficients of operators with ^Si quantum numbers, which 
yields a large correction in order v"^, yields a smaller correction in order v'^. 
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I. INTRODUCTION 



A formalism for the first-principles calculation of heavy-quarkonium decay rates in quan- 
tum chromodynamics (QCD) has been given in Ref. [1]. This formahsm is based on the 
effective field theory nonrelativistic quantum chromodynamics (NRQCD). In it, one can 
write the decay rate of a quarkonium state H as 



where F„ is a perturbatively calculable short-distance coefficient, m is the heavy-quark mass, 
the On are four-fermion NRQCD operators, and dn is the mass dimension of On- The terms 
in the sum over n may be classified according to their orders in v (Ref. [1]), where v is 
the heavy-quark-antiquark relative velocity. For charmonium, v'^ ~ 0.3; for bottomonium, 
^ 0.1. 

We concern ourselves in this paper with the decay of ^Sq quarkonium into light hadrons 
and the decays of ^Si quarkonium into lepton pairs and into light hadrons. The coefficients 
of the operators of leading order in v and of relative order have been computed previously 
[1, 2, 3, 4, 5, 6, 7, 8, 9]. Some of the coefficients of the order-f ^ operators are sufficiently large 
as to cast doubt on the convergence of the v expansion for charmonium and bottomonium. 
In particular, the order-v^ correction to the rate for the decay of ^5"! quarkonium into light 
hadrons is — 5.32(i>^), where {v'^) is the ratio of the expectation values of the order-v'^ and 
order-i>° operators in the quarkonium state. Hence, in the case of charmonium, the order-i>^ 
correction is more than 100%. 

In this paper, we compute the short-distance coefficients of the decay operators, through 
order f that have the same quantum numbers as the quarkonium. Our calculations confirm 
previous results for the short-distance coefficients of the order-v^ operators. We find that 
the V expansion is well behaved for the decays of ^Sq quarkonium and for the decay of 
^Si quarkonium into lepton pairs. In the case of the decay of ^Si quarkonium into light 
hadrons, large coefficients are associated with some of the operators of higher order in v. 
For the operators with "^S*! quantum numbers, a large correction to the decay rate appears 
in order v"^, but the correction in order is considerably smaller. This suggests that the 
V expansion for operators with a given quantum number may converge well once one goes 
beyond the first nontrivial order. 
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II. NRQCD DECAY RATES 



In this section, we present the NRQCD factorization expressions for the rates of ^5*0 
quarkonium {e.g. rjc or 775) decay to hght hadrons (LH), ^Si quarkonium {e.g. J/ip or T) 
decay to hght hadrons, ^Sq quarkonium decay to two photons, and ^Si quarkonium decay 
to e~^e~ . 

Through relative order v^, the rate for the decay of a ^5*0 state into hght hadrons is given 

by 



rCSo ^ LH) = 



^^^CSo\0,CSo)\'So) + ^^CS,\V,CSo)\'S, 



) 



+ ^^CSo\OseS,)\'So) + ^^CSo\OsCSo)\'So) 
+ ^^^CSo\OsCPi)\'So) + ^^^eSo\QlCSo)\'So) 

The operators appearing in Eq. (2.1) are defined by 

ViCSo) = 

^t(-iD)T,x-xH-tD)T,^, 
^t(_iD)2xx^(-iD)2^, 
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^W(-fD)2^ + ^t(_|D)2^^t^ 



Os{'S,) 
OsCSo) 
OsCPi) 
QlCSo) 



^^xx^D • + • D)^ -ij\B-gE + gE- D)xxV 
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(2.2c) 
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where the subscript 1 or 8 indicates that the operator is a color singlet or a color octet, 
the superscript labels the three dimension- 10 operators, ip is the Pauli-spinor field that 
annihilates a heavy quark, x^ is the Pauli-spinor field that annihilates a heavy antiquark, 
jjti — Qn igA^ is the gauge-covariant derivative, A is the SU(3)-matrix- valued gauge 
field, g is the QCD couphng constant, = G'°^ where = 9^^^^ - d"" Ai" + ig[Ai',A''] 
is the gluon field strength, and the cr' are Pauh matrices. The operator D is defined by 
X^Tiil^ — x^(D'0) — (Dx)^. The relative signs of the terms in each of these operators 



(and, in particular, ) are fixed by the requirements of hermiticity and charge-conjugation 
invariance. 

The matrix element of Qf does not appear in Eq. (2.1) because, as we show in Appendix A, 
it can be eliminated in favor of Q\ and Qf through the use of the equations of motion. From 
the velocity scaling rules in Ref. [1], we find that, in the ^5*0 state, the operator Ci(^S'o) has 
a matrix element of relative order v^, the operator Pi (^5*0) has a matrix element of relative 
order v"^, the operator Os{^Si) has a matrix element of relative order v^, and the operators 
Osi^So), OsC-Pi), Qli^So), Qli^So), and QI^Sq) have matrix elements of relative order v^. 

The contributions of order and order to the short-distance coefficient Fi(^S'o) have 
been computed in Refs. [2, 3] and are given in Ref. [1]: 



FiCSo) 



TlCp 2 



a;{2m){ 1 + 



\ ^ /199 1 37r2 



Ca - ^Uf 



a 



^ , (2.3) 



TT 



where A^^ = 3 is the number of colors, Cp = {N^ - l)/(2A'c) = 4/3, and Ca = N^. The 
contribution of order to Gi(^S'o) has been computed in Refs. [1,5].^ It is 



1 AttCfo 
Ci{ ^0) = rr^o; 



3A^c 



(2.4) 



We note that, to leading order in a^, [Gi(^So)]/[Fi(^So)] = —4/3. Hence, the first relativistic 
correction is sizable in the case of the rjc- The contributions of order and order to 
-^8('^'S'i) and Fs{^Sq) have been computed by Petrelli et al. [10]: 
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(2.5b) 



where is the QCD renormalization scale, n/ is the number of light-quark fiavors. Bp — 
(N^ - 4)/(4A^e) = 5/12, Tp = 1/2, and bo = (11/6)Ca - (2/3)Tpn/. The contribution of 



^ Short-distance coefficients can be extracted from the results in Ref. [5] by first making the substitution 
l/M^j,gQjj (l/4m^)(l — e/m), where — e is the binding energy, and then making the identification 
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order q;^ to Fs{^Pi) can be deduced from the results in Appendix A 2 of Ref. [1]: 



„ ,1 „ , TriVc 9 
FsCPi) = -^al. 



(2.6) 



Owing to energy conservation, the operators associated with the short-distance coef- 
ficients HK^Sq) and Hf(^So) cannot be distinguished from each other in the Born-level 
decay of on-shell quarks. Consequently, if one uses on-shell matching between NRQCD 
and full QCD to compute the short-distance coefficients in Born-level decay processes, one 
can compute only HK^Sq) + Hf(^So), not the individual coefficients. It is the quantity 
HlCSo) + Hf(^So) that we compute in this paper. 

Through relative order v"^, the decay rate for a ^5"! state into light hadrons is 



r{^Si LH) 



Gi(^-Si) 



+ ^^^{'S,\OsCSo)\'S,) + ^^^eS,\OseS^)\'S^) 



+ E 



Fsi^Pj) ,3q irn f3D \13q \ , -^l(^'^l) /3 c lr»lr3c M3( 



J=0,l,2 
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S,\Qi{'S^)\'S,). 



(2.7) 



The operator Os(^Sq) is defined in Eq. (2. 2d), and the operator O^i^Si) is defined in 
Eq. (2.2c). The remaining operators in Eq. (2.7) are defined by 



Oif^i) 
Pi('^i) 

0,{'Po) 

OsCPi) 

OsCP2) 



^ [z^Vx ■ x^(t(-|D) V + z/;to-(-|D)2x • xW 
^V^(-|D • tr) W(-iD • o-)7;V', 
i^t(_|D X cT)TaX ■ X^(-f D X o-)T,^, 
^t(_iB(v^))T„xx^(-iB(V^))T„V', 

V'M-|D)\-xM-iD)^v, 
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^to-x ■ xV(-|D)V + ^V(-|D)'ix ■ 
V'Vxx^ • (TiB -gE + gE- D)V - V^<^(D -gE + gE- D)x • X^o'V' 



(2.8a) 
(2.8b) 

(2.8c) 

(2.8d) 
(2.8e) 
(2.8f) 
(2.8g) 

(2.8h) 



The operator QK^Si) does not appear in Eq. (2.7) because, as we show in Appendix A, 
it can be ehminated in favor of Q\{^Si) and QK^Si) through the use of the equations of 



motion. Prom the velocity-scaling rules in Ref. [1], we find that, in the ^Si state, the operator 
Oi{^Si) has a matrix element of relative order v^, the operator Vi{^Si) has a matrix element 
of relative order v^, the operator Os(^So) has a matrix element of relative order v^, and the 
operators Os{'S,), OgCi^o), (^s^A), 08('i^2), Q\{'Si), Q?(^5i), and QH^S,) have matrix 
elements of relative order v^. 

The order-a^ and order-a^ contributions to the short-distance coefficient Fi{^Si) were 
computed by Mackenzie and Lepage [4] and can be found in Ref. [1], as can the order-a^ 
contribution: 
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where Q is the electric charge of the heavy quark, and the Qi are the electric charges of the 
light quarks. The order-a^ contribution to the short- distance coefficient Gi{^Si) is computed 
in Ref. [5]: 



5(197r^ - 132) 



(2.10) 
-5.32. Hence, 
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To leading order in al Gi{^Si)/[m^Fi{^Si)] = -{IQn^ - 132)/[12(7r2 - 9)] 
the relativistic correction to J/ ip decay is greater in magnitude than the leading contribution. 
This situation casts some doubt on the validity of the v expansion. Wc investigate this issue 
further in this paper by calculating corrections of relative order v^. The order-a^ and 
order-dig contributions to the short-distance coefficients Fs(^So) and Fs{^Si) are given in 
Eqs. (2.5b) and (2.5a), respectively. The order-a^ and order-a^ contributions to the short- 
distance coefficients F8{^Pj) have been computed by PetreUi et al. [10]: 
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(2.11c) 



where //a is the NRQCD renormahzation scale. The contribution to Fs{^Pi) of order van- 
ishes because Yang's theorem [11] forbids the decay of a spin-one particle into two equivalent 
massless vector particles (gluons). The contributions from decay into a light quark- ant iquark 
pair vanish because the ^Pj states are even under charge conjugation. Again, the individual 
quantities HK^Si) and Hf{^Si) cannot be distinguished in processes in which the heavy 
quark and antiquark decay on shell. We compute the quantity Hl{^Si) + Hf{^Si) in this 
paper. 

Through relative-order v^, the decay of a ^-So state into two photons is given by 



rCSo ^ 77) 



(oixvr-^o) 



+ 



Re 







+ ^^^^;^(^5o|^n-fD)^x|0)(0|x^(-|D)V|%) 



+ 



m 
^77 
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^0) 



Re 



C5o|V'^x|0)(0|x^(-iD)Vr^o) 



The product of matrix elements Re 



(2.12) 
, which is of 



(^'^olV'^xlO) (0|x^(D -gB + gB- T>)^\'So 
relative order v^, does not appear in Eq. (2.12) because, as we show in Appendix A, it can 
be eliminated in favor of the products of matrix elements in the last two terms of Eq. (2.12) 
through the use of the equations of motion. From the velocity-scaling rules in Ref. [1], we 
find that, in Eq. (2.12), the product of matrix elements in the first line is of relative order 
v'^, the product of matrix elements in the second line is of relative order v"^, and the products 
of matrix elements in the third and fourth lines are of relative order v^. 

The order-CK^ and order-a^ci;^ contributions to the short-distance coefficient Fj^(^So) are 
calculated in Refs. [2, 3, 6] and are given in Ref. [1]: 

^2 



TT 



1 + — - 5 Cf — 



a. 



TT 



(2.13) 



The order-CK^ contribution to Gjj{^So) is computed in Refs. [1, 5]: 



'^iiK '-'oj — 5 — . 



(2.14) 



To leading order in a^, G^^{} Sq) /[m^ F-y^{} Sq)] = —4/3. Hence, the first relativistic correc- 
tion to this process is substantial for the rjc- In this paper, we compute the combination of 
short-distance coefficients J^Sq) + J^Sq). 
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Through relative order v^, the rate for a ^Si state to decay into an e+e pair is 



r{^Si ^ e+e") 



-^ee(^'S'i) 



+ 



Re 



+ 



H. 



rrr 

1 (Z 



rrv 



S,|VV(-iD)\|0) . (0|xV(-iD)VI'5i) 



f5i|V^Vx|0)-(0|xV(-|D)V|'^i) . (2.15) 



The product of matrix elements Re 



which is of relative order f'^, does not appear in Eq. (2.15) because, as we show in Ap- 
pendix A, it can be eliminated in favor of the products of matrix clement in the last two 
terms of Eq. (2.15) through the use of the equations of motion. In Eq. (2.15), the product 
of matrix elements in the first line is of relative order the product of matrix elements in 
the second line is of relative order v^, and the products of matrix elements in the third and 
fourth lines are of relative order v^. 

The order-a^ and order-a^a^ contributions to the short-distance coefficient Feei^Si) are 
calculated in Refs. [7, 8] and are given in Ref . [1] . The order-a^a^ contribution is calculated 
in Ref. [9]. Altogether, these contributions give 



l-ACp 



as{m) 



TT 
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The order-a^ contribution to G^^Q-Sq) is computed in Refs. [1, 5]: 
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-a 



(2.16) 



(2.17) 



To leading order in ccg, Gee{^Si)/[m'^Fee{^Si)] — —4/3. Hence, the first relativistic correction 
to this process is substantial for the J/ip. In this paper, we compute the combination of 
short-distance coefficients Hl^i^Si) + Hl^i^Si). 



III. SPIN PROJECTORS 



In computing the quarkonium decay rates, we use the covariant spin-projector method 
[12, 13] to identify spin-singlet and spin-triplet amplitudes. For purposes of the computations 
in this paper, we need projection operators accurate at least through relative order v^. In 
this section, we compute the required projectors to all orders in v. 
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The Dirac spinors, with the standard nonrelativistic normahzation, may be written as 

„(p) = ( J_ ], (3.1a) 




"(-P) = ..+™ ■ I (3.1b) 



where ^ and 77 are two-component Pauh spinors, and E(p) — y^m? + p^. We take the heavy 
quark and antiquark momenta to be 

PQ = (l/2)P + p, (3.2a) 

Pq = (l/2)P-p, (3.2b) 

respectively, where in the quarkonium rest frame, 

P = {2E{p),0), (3.3a) 
p = (0,p). (3.3b) 

Using Eq. (3.1), it is straightforward to express spin-singlet and spin-triplet combinations 
of spinor bilinears in terms of Dirac matrices and to write them in a covariant form. In the 
spin-singlet case, we have 

no(P,p) = - J]ii(p,Ai)^;(-p,A2)(|Ai|A2|00) 

Ai,A2 

1 E + m f Q:-p\l + 7o / oc ■ p \ 

1 + 7^ n 75 1 - 7^ 7o 



■s/2 2E V E + m J 2 ' \ E + m 

(i f + m+ ^) {\f>-m-^), (3.4) 



2V2E{E + m) ' 4.E 

where ai and 7^^ are Dirac matrices in the Dirac representation, 75 = i7°7^7^7^, and we 
have chosen the normahzation so that the projector (3.4) corresponds in NRQCD to the 
projector I/V2, where / is a unit Pauh matrix.^ We note that E{p) may be written in a 
Lorentz invariant fashion as 

E{p) = (1/2)V7^. (3.5) 



^ In Eq. (3.4), the standard Clebsch-Gordan coefficients are appropriate if the spinors in Eq. (3.1) are 
related to each other through a unitary transformation, which preserves the SU(2) algebra, such as the 

charge-conjugation transformation rj = — i(T2^- One such choice of spinors ^^^~(gj (l) ^^'^ 



In the case of a spin-triplet state with polarization e, we have 
MP,p,e) = 5^tz(p,Ai)^;(-p,A2)(iAi|A2|le) 

Ai,A2 

1 E + mf^^a^\l+jo^_J^_^ ^^ 



^ 2E \ E + mJ 2 V E + m 

—1 P ^ IF 

(i /> + ^+^)Z:_L_ ^(i p-m-^). (3.6) 



2^/2E{E + m) AE 
Here, 1 1 e) is the rotationally invariant linear combination 1 1 e) = 1 1 1) — e+ 1 1 — 1) — €3 1 1 0) , 
with = (l/\/2)(ei ± 262). Wc have chosen the normalization so that the projector (3.6) 
corresponds in NRQCD to the projector cr • e/\/2. The expressions (3.4) and (3.6) are valid 
to all orders in v. 



IV. RELATIVISTIC CORRECTIONS TO ^Sq DECAYS 

In this section we compute the short-distance coefficients that appear in the corrections 
through relative order to ^5*0 quarkonium decays into two photons and into light hadrons 
(two gluons). 

We begin with the case of decay into two photons. We take the definitions of the heavy 
quark and antiquark momenta given in Eq. (3.2) and work in the quarkonium rest frame, as 
defined in Eq. (3.3). We take the outgoing photon momenta to be k and q, with polarization 
indices /i and u, respectively. Consider first the diagram in which the quark emits the photon 
with momentum k. The spin-singlet amplitude corresponding to this diagram is 

- + m 



A(sing^77) = -ie^Q^Tr 



1 



2y/2E{E + m) 2pQ ■ k 
2s/2E{E + m) 2pQ ■ k 



Tr 



-7 

-2pQ • k 



Tr 



Y ^1^{E + m+ f)^-^-t^{-E - m- f) 



(4.1) 



T] = 

convention is ^ 



for Ai 



or and 77= I or 

With this convention, the Clebsch-Gordan coefBcients in Eq. (3.4) must be multiphed by an additional 
factor (-l)(-i/2+A2). 



±1/2 and A2 
and T] ■ 



±1/2, respectively. On the other hand, a popular 
for Ai = ±1/2 and A2 = ±1/2, respectively. 
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where e is the electromagnetic couphng constant.^ In the projector (1 + 7o)/2 in the last 
line, the term proportional to 1 gives a vanishing trace, while the term proportional to 70 
gives 

771 1 

Ai(sing ^ 77) = (4-2) 
Similarly, the diagram in which the antiquark emits the photon with momentum k yields an 
amplitude 

771 1 

A,(sing ^ 77) = e'Q'-^e^^-^'k, -. (4.3) 

Adding Ai(sing — > 77) and A2(sing 77), we obtain the complete amplitude for ^Sq 
charmonium decay into two photons: 

We project out the 5"- wave part of the amplitude by averaging over the angles of p: 



^/2E " - |p||k|cos^ + |p||k| cos^ 



-e^Q^e^P^'^kp ^ In (4.5) 





IpI 


E- 


IpI 



where we have used = E^. Multiplying the expression (4.5) by its complex conjugate, by 
the two-body phase space l/(87r), and by a factor 1/2! for two identical particles in the final 
state, we obtain the decay width for a ^5*0 QQ state into two photons: 

r( ^0-77) = -^^ In ^3^. (4.6) 

Here, and in succeeding computations of the decay widths of two-particle states, we suppress 
a factor of the inverse volume that is associated with the normalization of the initial state. 

From Eq. (2.12), we find that the decay width for a ^5*0 QQ state into two photons in 
NRQCD in order a° and through relative order v'^ is 

rNRQCD('^o ^ 77) = 2[il/m^)F,,CSo) + (pVm^)G^-yC'5o) + {pVm')HU'So) 



77 ^ 

+ {pym')H'^^CSo)], (4.7) 



where the factor two on the right side of Eq. (4.7) comes from the spin factor for normalized 
heavy-quark states. 



^ In computing the short-distance coefBcients for the electromagnetic decay processes, we suppress trival 
color factors, which ultimately cancel when one matches decay rates in full QCD and NRQCD. 
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Comparing powers of p^/m^ in Eqs. (4.6) and (4.7), we obtain the short-distance coeffi- 
cients at leading order in as'- 

F^^CSo) = 27^Q^Q;^ (4.8a) 
G,,CSo) = -yQ'«^ (4.8b) 

H'^^CSo) + H',,CSo) = ^QW. (4.8c) 



Our results for ^^^(^5*0) and G^^(^So) confirm those given in Refs. [1, 2, 3, 6] and [1], 
respectively. Our result for H}^^{^So) + H^^{^So) is new. 

At leading order in ccs, the decay of a ^Sq QQ state to hght hadrons proceeds through 
an annihilation into two gluons. Hence, we may obtain the decay width for a ^5*0 QQ state 
into light hadrons by multiplying the width into two photons [Eq. (4.6)] by a color factor 
Cf/2 times al/{a^Q^): 

Prom Eq. (2.1), we find that the decay width for a ^Sq QQ state into two photons in NRQCD 
in order and through relative order v'^ is 

rNRQCDC-5o ^ LH) = 2N,[{l/m^)F,CSo) + {pym'')G,CSo) + {pVm'')HlCSo) 

+ m')HlC So)], (4.10) 

where the factor 2Nc on the right side of Eq. (4.10) comes from the spin and color factors 
for normalized heavy-quark states. The matrix elements of the color-octet operators do 
not contribute to Eq. (4.10) in order a'j.. Comparing Eqs. (4.9) and (4.10), we obtain the 
short-distance coefficients at leading order in as'- 

FiCSo) = (4.11a) 
G.CSo) = (4.11b) 
HlCSo) + HU^So) = (4.11c) 

Our result for Fi{^So) is in agreement with that given in Refs. [1, 2, 3], and our result for 
GiC'So) is in agreement with that given in Ref. [1]. Our result for Hl(^So) + Hf(^So) is new. 
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V. RELATIVISTIC CORRECTIONS TO ^Si DECAY TO e+e' 



Next we turn to the case of the decay of a ^Si quarkonium state into an e'^e~ pair. Again, 
we work in the quarkonium rest frame defined in Eq. (3.3). The amphtude for a quark and 
antiquark in a spin-triplet state with the momenta given in Eq. (3.2) to decay into a virtual 
photon with polarization index /j, is given by 



A(trip^7*) = ieQTV[ni(P,p,e)7^] 



Q\/2 



_E{E + m) 

We can project out the S-wsve part of the amplitude by averaging over the angles of p: 



(5.1) 



1 f 

A{^Si^-f*) = - / d{cos0)ieQV2 

2 J-i 



E{E + m) ^ 



(5.2) 



In order to obtain the decay rate into an e'^e" pair, we multiply the expression (5.2) by its 
complex conjugate with index u, by a photon-propagator factor —ignp/k^, by a complex- 
conjugated photon-propagator factor igua/k'^, and by twice the imaginary part of the e'^e~- 
pair contribution to the photon's vacuum polarization, namely, {gp^k^ — kpko^){—2/3)a. Here 
k is the virtual photon's momentum. The result is 



' ^ = ^E^[-3 + 3E 



(5.3) 



where we have used k ■ e — 0, e ■ e* — —1, and k"^ — 4£'^. 

From Eq. (2.15), we see that, in NRQCD through relative order v^, the decay width for 
a ^5*1 QQ state into an e+e" pair is 



/3 



NRQCD (' 5*1 



e e 



+ {pVm')HlCS,)]. (5.4) 



The factor two on the right side of Eq. (5.4) comes from the spin factor for normalized 
heavy-quark states. 

Comparing powers of p^ /rn? in Eqs. (5.3) and (5.4), we obtain the short-distance coeffi- 
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cients at leading order in ag'. 




(5.5b) 



(5.5a) 



(5.5c) 



Our result for Fee(^S'i) agrees with that given in Refs. [1, 7, 8], and our result for Ged^Si) 
agrees with that given in Refs. [1]. Our result for Hl^i^Si) + Hl^i^Si) is new. 

VI. RELATIVISTIC CORRECTIONS TO DECAY TO LIGHT HADRONS 

In the decay of a heavy-quark-antiquark state, diagrams in which only two of the final- 
state gluons attach to the heavy-quark line have a common heavy-quark color factor. Hence, 
(Abelian) charge-conjugation symmetry forbids such diagrams in the decay of a ^5"! state. 
Furthermore, color conservation forbids diagrams in which only one of the final-state gluons 
attaches to the heavy-quark line. Thus, in leading order in a^, a ^5*1 heavy-quark-antiquark 
state decays into three gluons, and the decay proceeds through diagrams in which all three 
gluons attach to the heavy-quark line. (Since no triple-gluon vertices appear, there are no 
ghost contributions.) 

In this decay process, in contrast with the decay processes that we have analyzed in the 
preceding sections, the kinematics allow one of the final-state gluons to have zero energy. 
Hence, the possibility arises that the decay rate contains an infrared (IR) divergence. Simple 
power counting arguments show that an IR divergence can arise only if the soft gluon 
attaches to an incoming (on-shell) heavy-quark or heavy-antiquark leg. Therefore, one can 
use NRQCD to analyze the interaction of this soft gluon with the heavy quark. 

One can see from power-counting arguments, that, through relative order t;^, in the 
Coulomb gauge, a gluon that interacts with a quark or an antiquark can yield an IR di- 
vergence only if the interactions arc of the type ip'^'D ■ Aip or x''"D ■ A^. (The ip'^'B ■ cri/j 
and x"l"B ■ crx interactions have the correct dimensions to produce an IR divergence, but the 
factors of B bring in powers of the gluon momentum that protect against an IR divergence.) 
The factor of D translates into a factor of the incoming quark or antiquark momentum. 
Factors of the gluon momentum do not appear since they are orthogonal to the gluon prop- 
agator in the Coulomb gauge. Therefore, the interactions of the gluon yield two factors of 
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the incoming quark or antiquark momentum in the squared amphtude. Two additional fac- 
tors of the incoming quark or antiquark momentum are required in order to have a nonzero 
overlap with an incoming 5'-wave state. Hence, an IR divergence in the decay rate must be 
associated with at least four factors of the incoming quark or antiquark momentum. That 
is, an IR divergence can first appear in relative order v^. Because the soft gluon in a D • A 
interaction changes the incoming 5'-wave color-singlet quark-antiquark state into a P-wave 
color-octet quark-antiquark state, we expect that the IR divergence will be absorbed into 
matrix elements of the P-wave color-octet operators in Eq. (2.7). 

Now let us turn to the actual computation of the rate for a ^5"! QQ state to decay into 
three gluons. We present only the outlines of that calculation here. We used the symbolic 
manipulation program MATHEMATICA and the package FEYNCALC [14] to handle the tedious, 
but straightforward, details of the algebra. 

We regulate the anticipated IR divergence by computing in D = 4 — 2e dimensions. We 
work in the quarkonium rest frame, assign the incoming quark and antiquark momenta 
as in Eqs. (3.2) and (3.3), and take the outgoing gluon momenta to be ki, k^, and k^. 
First we compute the sum of the six Feynman diagrams for this process, making use of the 
projector (3.6). Although we are working m D = A — 2e dimensions, we can follow the 
approach of Ref. [10] and simply use the D-dimensional version of the spin-1 projector (3.6). 
As explained in Ref. [10], we need not consider projectors for the higher-spin evanescent 
NRQCD operators that appear in D dimensions because the contributions that contain an 
IR pole in one loop do not mix the higher-spin operators with the spin-1 operators. 

At this point, we could square the amplitude, integrate over the phase space, and expand 
in powers of p in order to obtain the desired result. However, the amount of algebra would 
be greatly reduced if we could make the expansion in powers of p before carrying out the 
phase-space integration. Such a strategy is complicated by the fact that the phase space 
depends on p through the total energy of the incoming QQ state, but we can make that 
dependence explicit by introducing a rescaling of phase-space integration variables: 

h ^ hE{p)/m. (6.1) 

Then, the final-state phase space transforms as 
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where 





{D-2f/D 




2 - 


m 




m 





5 E\p) 
1 - -e log 



(6.3) 



All of the dependence on p on the right side of Eq. (6.2) is contained in the explicit factor 
fip). The remaining factors correspond to the phase space evaluated at the QQ threshold 
p — 0. Therefore, after rescaling the ki according to Eq. (6.1), we can obtain the necessary 
expansion in powers of p by expanding the amplitude, its complex conjugate, and f(p) in 
powers of p before carrying out the phase-space integration. Note that an IR pole in e first 
appears in the rate, excluding the factor f{p), only in the relative-order v^. Hence, we can 
drop the term proportional to e on the right side of Eq. (6.3), which contributes an additional 
factor v'^. 

We expand the amplitude in a power series in p, through order p^. The terms containing 
no powers of p yield a pure S-wecve contribution. For the terms containing two powers of p, 
we extract the S-wave contribution by making the replacement 



PfiPu 



(6.4) 



For the terms containing four powers of p, we extract the S'-wave contribution by making 
the replacement 



Here, 



PiiPuPpPa P^T^^pa- 



T = n 



D-l 



T, 



and 



nupa — _ i^^jj _|_ [n/i;/Hp(^ + H^pHj,(^ + H^(jHj,p] , 

P P 



9p.v ~l~ 



(6.5) 
(6.6) 
(6.7) 
(6.8) 



4^2 (p) 

Next, we multiply the amplitude by its complex conjugate. We evaluate the gluon polariza- 
tion sums using the Feynman-gauge expression 



^p,^v Qp^f-i 

we evaluate the spin-triplet-state polarization sum using 



(6.9) 



(6.10) 
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and we divide by — 1 to obtain the average over the spin-triplet-state polarizations. Owing 
to the charge-conjugation invariance of the amplitude, only the part of the color factor that 
is symmetric in the color indices survives. It is given by 

1 MzilM^. 



Multiplying by this color factor and by f{p), we obtain the "squared matrix element" that 
must be integrated over the p — three-body phase space to obtain the decay rate. We 
write the coefficients of p°, p^, and p"^ in terms of the invariants 

s ^ 2ki- k2, (6.12a) 

t = 2ki- ks, (6.12b) 
u ^ 2k2- k3, (6.12c) 

where, since we have set p = in these coefficients, the energy-momentum conservation 
relation now reads 

ki + k2 + ks^2m. (6.13) 

The expressions for these coefficients in four dimensions are given in Appendix B. 
We re- write the coefficients of p°, p^, and p^ in terms of the invariants 



Xi 



2P ■ ki 



{2my 
It follows that 



(6.14) 



p=0 



s — 2m'^{xi + X2 — X3), (6.15a) 
t ^ 2m'^{xi- X2 + X3), (6.15b) 
u — 2m^(— Xi -|- X2 + X3). (6.15c) 

The £)-dimensional three-body phase space for decay of a particle of mass M is [10] 



'2c 3 3 

(S) ^rx^n(i--0-^^-^'^(2-E-0- (6-16) 



"^'^ 2(47r)3 

The phase space at p = is obtained by making the identification M — 2m. 
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It is convenient to make a further change of variables, so that the hmits of integration 
are independent of the integration variables. To this end, we write 



xi — X, (6.17a) 
X2 = i — xy, (6.17b) 
x^ = 1 — x{l — y). (6.17c) 

This change of variables is particularly useful in analyzing the infrared singularities, since 
it avoids the difficulty that, at the singular points Xi = 0, the range of integration in one of 
the variables Xj {i ^ j) vanishes. Now the phase space is 

where x and y range from to 1. 

In the cases of the terms proportional to p° and p^, the integrations over the phase space 
are IR finite, and we can carry out the integrations with D — A. Multiplying by 1/3! for 
three identical particles in the final state, we obtain 



In the case of the term proportional to p^, we must first separate the IR singular parts 
in the matrix element squared.^ These are given by 

= £;M*i) « 1(3 - - , - 2(2 - e),(l - 

1 1 1 



x"^ (1 — xy)"^ [1 — x{l — y)]' 



(6.21) 



In Ref. [10], an alternative method for dealing with the singular part was employed. The region of 
integration was partitioned into three regions that are related by interchange of the three gluon momenta. 
Only the region containing the singularity at a; = was retained, and the contribution from this region 

was multiplied by three to obtain the complete result. Th(^ method that wc present in this paper has the 
advantage that the limits of integration arc simpler, and lionco, the integrals arc evaluated more easily. 
Also, certain terms that cancel between the singular and non-singular contributions in the method of 
Ref. [10] never appear in the present method. 
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Integrating Mir over the phase space (6.18) and multiplying by 1/3! for three identical 
particles in the final state, we obtain 



r(4) _ 

^ IR — 



iN!-im-^) ,l /4vry- r^(l-6) (1-6)^(7-46) 
'e J ^ T\2-2e) {3 - 2e)^ 



m 



(6.22a) 



Neglecting terms of order e, we may write this expression as 



.(4) 
IR 



(iVe'-l)(^^c'-4) 



-a. 



7 /47r 
81e \ JP 



2e 



44 



6, (l-67^)r(l-6) _ 

^ r(2 - 26) 243 



(6.22b) 



where 7^ is Euler's constant. After subtracting the IR singular terms (6.21) from the 
integrand, we can carry out the phase-space integration over the remainder with D — A. 
Multiplying by 1/3! for three identical particles in the final state, we obtain 



p(4) _ (iVe^-l)(Ar,2-4) 3 
^ finite _fi Ar9 s 



rrr 



3563 1609 2 



2430 6480 



(6.23) 



The complete decay width in full QCD of a ^Si QQ state into fight hadrons through order 
v"^ is then 

^es^ ^ LH) = r(°) + r(2) + + r^tL, (6.24) 

where the quantities on the right side are given in Eqs. (6.19), (6.20), (6.22), and (6.23). 

To determine the short-distance coefficients, we match these results with the NRQCD 
expression for the decay width (2.7), evaluated in the ^Si QQ state. Since we have computed 
the full QCD decay rate in order a^, we must evaluate each contributing term in Eq. (2.7) 

with accuracy a^. The coefficients Fs(^Sq), Fs{^'Si), and Fs{^Pj) are of order Therefore, 
we must evaluate the corresponding matrix elements through order ag- We evaluate the 
matrix elements corresponding to the unknown coefficients Fi{^Si), Gi{^Si), Hl{^Si), and 
Hf{^Si) at order q;°. 

The color-octet matrix elements {^Si\OsCSo)\^Si) and (^5'i|C8(^5'i)p5'i) have a vanishing 
contribution at order q;° in the color-singlet QQ state. The order-aj contribution comes 
from four diagrams in which a gluon connects an initial-state Q oi Q with a final-state Q 
or Q. The interaction of the gluon with the Q or Q cannot be of the p ■ A form, since 
that interaction changes the orbital angular momentum by one unit. Any other NRQCD 
interaction must involve at least one power of the gluon momentum. Hence, it is easy to 
see, by simple power counting arguments, that the integration over the gluon momentum is 
ultraviolet (UV) power divergent. It therefore vanishes in dimensional regularization. 
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The color-octet matrix elements {^Si\Os{^Pj)\^Si) also have a vanishing contribution at 
order in the color-singlet QQ state. Again, the order-aj contribution comes from four 
diagrams in which a gluon connects an initial-state Q or Q with a final-state Q or Q. The 
contribution at leading order in v arises from p ■ A interactions between the gluon and the 
Q or Q. A straightforward computation yields 



{'s,\osepj)\'s,) = 



p^ 8{2J+1)Cf 



-a.c 



dk 
T' 



(6.25) 



This integral has logarithmic IR and UV divergences. Since it is scale invariant, it vanishes 
in dimensional regularization. It can be written as 



8l7r 



,,2e 



a. 



euv 



+ log(47r) - 7s 



h log(47r) - 7e 



(6.26) 



where euv and eiR are (4 — D)/2, and /i^ is the NRQCD renormalization scale. We renormal- 
ize the expression (6.26) in the modified minimal subtraction (MS) scheme by subtracting 
the contribution proportional to l/euv + log(47r) — 7e. The renormalized matrix element is 



{'Si\Os{'Pj)\'S,] 



p4 A{2J+l)CFfi 



2e 



MS 



Sin 



+ log(47r) - 7e 



(6.27) 



where we have made the identification ei^ — e. 

Making use of these results for the matrix elements, we find that the decay width in 
NRQCD in order is 



rNRQCD('5i ^ LH) = 2N,[(l/m')F,{'S,) + (p' /m'')G,('S,) + {p' lm')H]^{' S,) 



ri A. 



+ (pVm6)iJ^,(i5o)] + (pVrn^) E cji^sC^^)' (6-28a) 

J=0,l,2 

where the factor 2A^c in front of the square brackets comes from the color and spin factors 
for normalized heavy-quark states and 



cj 



,2c 



SIttN, ^ Va 



1 



+ log(47r) - 7b 



(6.28b) 



The short-distance coefficients F8(^Pj) have been computed D = 4 — 2e dimensions in 
order by Petrelli et al. [10]: 



4e 



r(i 



1 



FsCP2) 



0, 



r(2-2e)3-2e' 

r(l-e) 6-13e + 4e2 
r(2-2e) (3-2e)(5-2e)' 



(6.29a) 
(6.29b) 
(6.29c) 
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It follows that 

^ cM( P,) - - „. - j 



J=0,l,2 



(1 - e7E)r(l - e) 1 
r(2 - 2e) 15 



(6.30) 

where we have neglected terms of order e. 

Using Eq. (6.30), we can compare the width in full QCD [Eq. (6.24)] with the width in 
NRQCD [Eq. (6.28)] to compute the short- distance coefficients. As expected, the IR poles 
in e cancel, and we obtain 

c.e.0^M^(n_|>..)„3 ,,,, 

, „2,3c> _ {N'l-mNl-t)( 833 , 1609 , , 7 , 2m\ , 



(6.31c) 

Our result for Fi(^S'i) agrees with that given in Ref. [4], and our result for Gi(^Si) agrees 
with that given in Ref. [5]. Our result for HK^Si) + HK^Si) is new. 



VII. DISCUSSION 



In this paper, we have computed short-distance coefficients for the decays of a Sq heavy- 
quarkonium state to two photons and to light hadrons and the decays of a ^5"! heavy- 
quarkonium state to a lepton pair and to light hadrons. Specifically, wc have computed 
the coefficients of the operators whose matrix elements are of order and whose quantum 
numbers are those of the quarkonium state. 

In our computation, we are able to obtain only the combinations H^i^^^^Lj) -\- 
if^(^'^+^Lj), rather than the individual coefficients if^(^'^+^Lj) and if^(^'^+^Lj), because 
the corresponding operators, Q^i^^^^Lj) and Q^(^'^+^Lj), have identical matrix elements 
for on-shcll heavy quarks in the ccntcr-of-momcntum frame. In order to obtain the vahics of 
the individual coefficients, it would be necessary to consider matrix elements of the operators 
and ^i^^^^Lj) in which the heavy QQ interact with additional quanta before 
reaching the annihilation vertex. Alternatively, one could consider matrix elements of the 
operators Q^i^^'^^Lj), which, as we have shown in Appendix A, are related to the operators 
Qi^25+i^^^ and Q^i^^'^^Lj) through the equations of motion. 
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In Tables I-IV, we show the numerical values of the short-distance coefficients that appear 
through order for the decays that we consider in this paper. For each coefficient, we take 
into account only the contribution that is leading in a^. In each case, we normalize the short- 
distance coefficients to the coefficient of the operator whose matrix element is of leading order 
in V. In the third column of each table, we use the velocity-scaling rules [1] to estimate the 
size of the matrix element of the operator that is associated with each coefficient, relative to 
the size of the matrix element of leading order in v. In the case of the color-octet operators, 
we adopt the approach of Ref. [10], multiplying the velocity-scahng estimate by a factor 
l/{2Nc) to account for the relative spin and color normalizations of the color-singlet and 
color-octet operators as we have defined them in this paper. 

Coefficient Value Matrix Element 

F^^(^5o) 1 1 

G^^i^So) -1.33 ^2 

H}^^{^So)+H^^{^So) 1.51 

TABLE I: Short-distance coefficients and estimates of sizes of corresponding matrix elements for 
the decay of a quarkonium state to two photons. 



Coefficient Value Matrix Element 
Fii^So) 1 1 

Gii^So) -1.33 ^2 

Fsi^Si) 0.75nf v^/{2N,) 
FsCSo) 1.88 vy{2N^) 

Fs^Pi) 1.13 vyi2Nc) 



HlCSo) + HfCSo) 1.51 

TABLE II: Short-distance coefficients and estimates of sizes of corresponding matrix elements for 
the decay of a quarkonium state to Ught hadrons. 

In the case of charmonium, v"^ ~ 0.3 and as{mc) ~ 0.35. Then, we see from Tables I-III 
that the convergence of the v expansion is reasonable for the ^Sq decays into two photons 
and into light hadrons and for the ^Si decay into hght hadrons. 

On the other hand, the coefficients in Table IV cast some doubt on the convergence of the 
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Coefficient Value Matrix Element 

Feei^Si) 1 1 

GeeCSi) -1.33 

Hl{^S,)+Hl{^S,) 1.61 

TABLE III: Sfiort-distance coefficients and estimates of sizes of corresponding matrix elements for 
the decay of a ^Si quarkonium state to a lepton pair. 



Coefficient Value Matrix Element 

Fii^Si) 1 1 

Gi(35i) -5.32 ^2 

F8('5o) 11.647rK ^V(2iVc) 

Fsi^Sx) 4.66n/7r/a, v^/{2N^) 

Fsi^Po) 34.937r/a, vy{2N,) 

Fsi^Pi) 2.26 - 6.90n/ + 5.17n/ log(2m///A) i;^/(2Ar^) 

F8(^J'2) 9.3l7r/a, vy{2N,) 

Hli^Si) + H^i^Si) 7.62 + 1.791og(2m//XA) 



TABLE IV: Short-distance coefficients and estimates of sizes of corresponding matrix elements for 
the decay of a ^5*1 quarkonium state to light hadrons. 

V expansion in the case of the ^Si decay into light hadrons. In the case of charmonium, all of 
the contributions of higher order in v are larger in magnitude than the order-f " contribution, 
with the exception of the Hl{^ Si) + Hf{^ Si) contribution. The color-octet coefficients, other 
than Fs{^Pi), are enhanced by n/ag, relative to since the corresponding color-octet 

Fock states can decay into two gluons or into light-quark pairs, rather than into three gluons. 
In addition to this enhancement, some of the coefficients of ir/ag are quite large. However, 
one can, through a redefinition of the color-singlet operators, incorporate the factors 1 / {2Nc), 
which we have associated with the matrix elements, into the short- distance coefficients [10]. 
Then, aside from the tt/qs enhancement, only Fs{^Po) is especially large. In the case of the 
color-singlet coefficients, Gi{^Si) is quite large in magnitude relative to Fi{^Si). However, 
the quantity Hl(^Si) + Hf{^Si) is not significantly larger in magnitude than Gi{^Si), giving 
some hope that the v expansion may ultimately be well behaved. 
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The estimates of the sizes of the relativistic corrections strongly suggest that, in order to 
carry out a meaningful phenomenological analysis of S-wave quarkonium decays into light 
hadrons, one would need to take into account contributions beyond leading order in v. (For 
a further discussion of this point, see Ref. [15].) All of the contributions listed in Table IV, 
except for that of F^i^Pi), would be needed to achieve a precision of better than 50%. 

Unfortunately, most of the required matrix elements are unknown. However, the number 
of unknown quantities can be reduced drastically by making use of the heavy-quark spin 
symmetry and the vacuum-saturation approximation [1] , although the accuracy of these ap- 
proximations is not always sufficient to allow a calculation of the decay rates through relative 
order v^. Owing to the heavy-quark spin symmetry, the matrix elements of Oi(^So), Vi{^So), 
Osi^Si), Osi^So), Q\i^So), and QK^Sq) in a ^5*0 state are equal to the matrix elements of 
Ci('^i), Pii^Si), OsCSo), Osi^Si), Qli^Si), and QH^Si) in a state, respectively, up to 
corrections of relative order v^. Also owing to the heavy-quark spin symmetry, the matrix 
elements of the operators Os{^Pj) in a ^Si state are equal (2 J -|- l)/9 times the matrix 
element of Os{^Pi) in a ^5*0 state, up to corrections of relative order v'^. According to the 
vacuum-saturation approximation, the matrix elements of the operators for the electromag- 
netic decays are equal to the matrix elements of the color-singlet hadronic-dccay operators 
with the same quantum numbers, up to corrections of relative order f ^. It also follows from 
the vacuum-saturation approximation that the matrix element of Q\{'^^~^^Sj) is equal to the 
square of the matrix element of Vi{^^'^^Sj) divided by the matrix element of Oi(^^'^^Sj), 
up to corrections of order v^. However, the matrix element of QK'^^'^^Sj) is not known to 
be related to the others. 

The matrix elements of Oi{^Si) in the J/ip and T states are known from phenomenol- 
ogy. The matrix elements of Oi{^Si) and Vi{^Si) in the J/ip and T states have also been 
computed on the lattice [16], although the lattice determinations of the matrix elements of 
Vi{^Si) are rather imprecise, owing to large uncertainties in the perturbation series that 
relates the lattice and continuum matrix elements. According to the Gremm-Kapustin re- 
lation [17], for dimensionally regulated matrix elements, the matrix element of Vi{^Si) is 
equal to the matrix element of Oi{^Si) times (M — 2mpoic)/'ni, up to corrections of relative 
order v^. Here, M is the quarkonium mass, and mpoie is the heavy-quark pole mass. The 
remaining unknown operator matrix elements could, in principle, be determined in lattice 
numerical simulations. 
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APPENDIX A: RELATION BETWEEN THE OPERATORS OF ORDER 

In this appendix, we demonstrate that the operators Q\(^Sq) [Eq. (2.2)] are related 
to each other by the equations of motion, as are the operators Q\{^S]) [Eq. (2.8)], the 
vacuum-saturated versions of the Q\ {}Sq), and the vacuum-saturated versions of the Q\ {^Si). 
Wc assume that these operators are integrated over all space-time, so that we can employ 
integration by parts in re-writing them. 

We begin by considering the operator 

QlQSo) = ^\-lT)fxx\-{T)fi^. (Al) 

Now, 

X^(-iD)V = x\-id-gP^ + id-gA){-id-gA + id-gA)i^ 
= X^[2{-id - gAf + 2{id - gAf - {id + idf]il: 
= X^[Am{id^ - gAo) + 4m(i9o - gA^) - {id + ia)> 
= [Aimdo - {-idfy^, (A2) 

where we have used the equations of motion at leading order in v in the third line. Further- 
more, under integration by parts, which is equivalent to energy-momentum conservation in 
momentum space, 

[ip^{-iByx][4:imdo - {-idf]x^ilj {[-4mido - {idf][ilj^{-iByx]}x^ip. (A3) 

Let us focus on the first QQ bilinear on the right of Eq. (A3). It is 

[-4mido - {idf][ilj\-iBfx] = ijj^{-4m{ido + gAo){-iT>f - 4:m{-iBf{ido - gAo) 

+ Am[gAo, (-iD)2] - Am[ido, (-iD)^] 

-{id + idy{-iBy}x. (A4) 
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Then, using the equations of motion, we have 



+ 4m[^A, {Sf\ - 4m[i9o, {Sf\ 

- {id + idf{-iY) f}x 
= V'H(-^D)^ + 4m[(-^Do),(-^D)2]}x 
= ^^[{-iT)Y - 8m(D • + • D)]x, (A5) 

where, in arriving at the second equahty, we have dropped some terms proportional to 
[Dj, Dj\ — —2ig€ijkBk that are order v relative to the terms that we have retained. 
Thus, taking into account both QQ bilinears, we have 



QlCSo) ^ Ai-I^Y - (W2)(D -gB + gB- D)]xxV- 



(A6) 



Carrying out this procedure symmetrically on the left and right QQ bihnears of Q\(^So), 
we conclude that, under the equations of motion and integration by parts. 



QiCSo)^QiCSo) + {m/2)QtCSo) 



(A7) 



A similar analysis in the spin-triplet case yields 



12 /S 



.3/3 I 



(A8) 



The vacuum-saturated versions of these relations, which are relevant to the electromagnetic 

decays are 



and 



Re 



Cso\i^h\o){o\x\-fimso) 



+ (m/2) Re (^5o|V'^x|0>(0|x^(D • ^E + ^E • D)V'r^o) 



(A9) 



{'S,\i;^a{-Sfx\0) ■ (0|xV(-|D)V|'^i) 



Re 



f5i|V'Vx|0)-(0|xM-|D)V|'^i) 



+ (m/2) Re f 5i|^/;Vx|0) ■ (0|xV(i3 • ^E + ^E • Bj^l^^i) . (AlO) 
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APPENDIX B: SQUARED MATRIX ELEMENTS FOR DECAY INTO 
LIGHT HADRONS 

In this appendix we give the expressions for the terms of order p(°) , p*^^) , and p^^^ in the 
square of the matrix element for a ^Si QQ state to decay into hght hadrons (three gluons). 
These terms are denoted by and Ai^^\ respectively. The quantity m is the 

heavy-quark mass. The invariants s, t, and u are defined in Eq. (6.12). 



^,„, ^ {N!-im-i)2mS^ai _ ^^,^3 ^g^,^^ _ ^2^,^,^ ^ ^^3^ 

3 

+ 16mV - 12m^st^ + Ssh"^ - Smh^ + 2st^ + t^) 



j^i2) ^ _p2 (^c 1) (iV, 4) 81927r a, ^^^3^10^2 _ ^^^^8^3 _ 43^6^4 ^ ^^^4^5 

9 

- 2w?s^ + 512m^°st - 128m^s^i - 2bQm^sH + l2Qm^sH - ISrri^sH + sH 

- 256mSi^ + 176mS^t^ - AAm^sH^ + 5sH^ - A8mH^ + 120m^st^ - 37m^sH^ 
+ BsH'^ + 2AmH^ - ISm^st^ + 3sH^ - 2mH^ + st^) 

^ m?{Am? - sf{Am? - tf{s + tf ' ^^^^ 
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+ 223232m^S^ + 47104m^s^ - 14080mS^ + 880m^s^ + 1310720m^Si 
+ 2588672m^S2i - 4034560m^^s^i + 1176576m^°s^i + 97536m^s^i 

- 86848771^5^^^ + 12352777^s^t - 552777^5^^ + 1966080777^''t2 + 2588672777^^5^^ 

- 3047424777^^5^*^ + 984064777^°s^i^ + 371456777^s^i^ - 275264777^s^t2 

+ 56032777S''t^ - 4688777^sV + 155sV + 1409024777^^t^ - 4034560777^^si^ 
+ 98mAm^°sH^ + 525312m^sH^ - 408960777^5^^^ + 108304771^5^^^ 

- 12452777^5 V + 620s V - 1445888777^^^^ + 1176576777^°si^ + 371456777^s^i^ 

- 408960777^=^^^ + 128640777 Vi^ - 18468777^s^i^ + 1240s¥ + 223232777^°^^ 
+ 97536771® st^ - 27b2Um^sH^ + 108304777^s2t^ - 1846877i2s^t5 ^ 155055^5 
+ 47104777®^^ - 86848777^5^^ + 56032777^^^^ - 12452777^5^^^ + 1240s¥ 

- 14080777^^^ + 12352777^^^ - 4688777^s^t^ + 620sV + 880777^^^ - 5527772si^ 

+ 15552^8)——— -J— — -. (B3) 

^ 777^(47772 - s)4(47772 - tY{s ^ ^ 
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